Finite range effects in (d,p) reactions 
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With the increasing interest in using (d,p) transfer reactions to extract structure and astrophysical 
information, it is important to evaluate the accuracy of common approximations in reaction theory. 
Starting from the zero-range adiabatic wave model, which takes into account deuteron breakup 
in the transfer process, we evaluate the importance of the finite range of the n — p interaction in 
calculating the adiabatic deuteron wave (as in Johnson and Tandy) as well as in evaluating the 
transfer amplitude. Our study covers a wide variety of targets as well as a large range of beam 
energies. Whereas at low beam energies finite range effects are small (below 10%), we find these 
effects to become important at intermediate energies (20 MeV/u) calling for an exact treatment of 
finite range in the analysis of (d,p) reactions measured at fragmentation facilities. 

PACS numbers: 24.10.Ht; 24.10.Eq; 25.55.Hp 



I. INTRODUCTION 

Since the early days of nuclear physics, transfer reac- 
tions have been a preferred tool to study spectroscopic 
properties of nuclei and have been widely used to de- 
termine single particle structure across the nuclear chart 
(e.g. [il-01). Such studies have allowed a better under- 
standing of detailed features of nuclear interactions. It 
was through a systematic study of single nucleon trans- 
fer on Sn isotopes, that we now understand the reduction 
of the spin-orbit strength when moving toward neutron 
rich systems [4]. One would also like to use the transfer 
reaction method to discriminate between effective inter- 
actions used in the shell model, as suggested in the study 
of (d,p) on Ni isotopes Q. Nowadays the spectroscopy 
of exotic nuclei can be studied using inverse kinematics. 
Pioneering studies |10l - [l4| hold promise for applying this 
technique more broadly, especially in new generation rare 
isotope facilities, where beam rates will be enhanced. 

Another intriguing aspect of nuclear structure is the 
role of correlations. The independent particle shell model 
of course neglects all correlations. State-of-the-art shell 
model include some correlations effectively through the 
residual interactions. Electron knockout experiments 
have shown a 30% decrease in the spectroscopic factors 
of closed shell nuclei as compared to the independent 
shell model predictions, a reduction that is understood in 
terms of short range correlations [la] . Reduction factors 
from nuclear knockout experiments can be much larger, 
depending on the difference between the neutron and 
proton separation energies [l7| . However, the physical 
reason for such a large suppression is still unclear, long 
range correlations being a possibility [T^]. What is most 
intriguing is that spectroscopic factors from transfer re- 
actions show no dependence on the difference in neutron- 
proton separation energies. Large surveys of ground state 
spectroscopic factors from (d,p) reactions p], including 
nuclei with a wide range of separation energies (0.5-19 
MeV) show good agreement with large scale shell model 



using modern effective interactions. 

The reconciliation of these results with those from 
knockout experiments [15j is proving difficult and is the 
subject of recent work on new approaches to the calcula- 
tion of overlap functions [la, [l^. The resolution of these 
important physics questions relies on the accuracy of the 
reaction model being used. It is thus of paramount im- 
portance that the reaction theory is well founded and 
uncertainties are well understood. In this respect there 
have been a number of works looking at specific aspects of 
the reactions theory used for the analysis of (d,p) transfer 
data (uncertainties in the optical potentials (24| , coupling 
to excited states of the target [2^ and ambiguities due to 
the single particle wave functio ns |26l. [27| and new ways 
of calculating overlap functions [ISl . Il9| ) . In this work we 
explore another aspect: the consequences of the non-zero 
range of the n—p interaction that plays a key role in the 
theory is several different ways. 

Over the past four decades of work on (d,p) reactions, 
different approximations were made, not only regard- 
ing the optical potentials, the deuteron and final single- 
particle wavefunctions, but also in the reaction mecha- 
nism and the evaluation of the transfer matrix element. 
It may thus be confusing to realize the large range of 
spectroscopic factors available in the literature. System- 
atic and consistent studies 0, Q-Q provide a much better 
overall assessment of the situation. In [^, 0J an extensive 
survey of ground state spectroscopic factors from charge 
Z = 3 — 24 was performed using the same reaction model, 
the same global potentials and the same single particle 
parameters. With these same assumptions, spectroscopic 
factors extracted from 235 sets of (d,p) data were found 
consistent with shell model predictions to within ±20%. 
An identical study was performed on excited states [8|, 
and although there were a few unresolved cases (such as 
the Ni isotopes), the overall agreement with shell model 
was of the order of 40%. Ni isotopes were studied sepa- 
rately [3| and a better overall description of these nuclei 
was found with a modified effective interaction in the 



shell-model calculations. 

All these studies rely on the adiabatic distorted- 
wave approximation (ADWA) developed by Johnson and 
Soper '28'] and the local energy approximation (LEA) |2<l| 
to take into account the finite range of the n — p interac- 
tion Vnp in evaluating the (d, p) transition matrix. 

In 28], a three-body theory for A(d,p)B was developed 
taking into account deuteron breakup which is known to 
be important even for reactions on stable nuclei. In [2^ 
the zero range approximation is made for the n—p inter- 
action. Then, the transfer matrix element reduces to a 
form similar to the zero range Distorted Wave Born Ap- 
proximation (DWBA) where the scattering wavefunction 
in the incident channel is calculated with an adiabatic po- 
tential consisting of the sum of the proton and neutron 
potentials evaluated at half the deuteron energy, instead 
of the deuteron optical potential extracted from (d,d) 
data. ADWA typically decreases the radius and diffuse- 
ness of the distorting potential |30| compared with typ- 
ical values for deuteron optical potentials deduced from 
elastic deuteron scattering and often provides a better 
description of the data [31.] . 

An adiabatic theory including finite range effects was 
formally developed by Johnson and Tandy [32|. Using 
a Weinberg expansion in the deuteron channel, Johnson 
and Tandy arrive at a set of coupled channel equations to 
describe the relative motion of the centre of mass of the 
neutron and proton relative to the target. The solution of 
the coupled channels equations gives a three-body wave- 
function that is a coherent superposition of the bound 
(deuteron) and break-up continuum states of the n—p 
system. 

In this paper we will confine ourselves to the simplest 
version of this theory in which only the first term in the 
Weinberg expansion is retained. This assumes that the 
only significant break-up components in the three-body 
wavefunction have sufficiently small energies that inside 
the range of Vnp, the relevant n—p scattering states 
have the same radial shape as the deuteron ground state 
wave function. When only the first Weinberg term is in- 
cluded, the coupled equations reduce to an optical-model- 
like equation for the three-body wavefunction where the 
distorting potential is the sum of the neutron and proton 
optical potentials multiplied by the neutron-proton inter- 
action, folded over the deuteron bound state. The effect 
of other components |33| have been shown to be signifi- 
cant at Ed = 88 MeV but their effects at lower energies 
are unknown. 

In the early sixties, it was already understood that fi- 
nite range effects were important in (d,p) reactions, how- 
ever due to computational limitations, all calculations 
were performed in zero range. A very popular proce- 
dure to correct a zero r ang e calculation was developed 
by Buttle and Goldfarb [231 ■ The standard implementa- 
tion of this method, the so-called local energy approxi- 
mation (LEA) , takes only the first term of the expansion 
presented in 29]. For deuteron energies well above the 
Coulomb barrier, it is not clear that this procedure is suf- 



ficiently accurate. A simple estimate of the modification 
of the Johnson and Soper potential due to a finite range 
Vnp were reported in |37| . 

In this work, we perform a systematic study of finite 
range effects in (d,p) reactions within the framework of 
ADWA. We consider 26 (d,p) reactions on stable targets, 
involving nuclei with masses ranging A = 12 — 208 and 
deuteron energies Ed = 2 — 70 MeV. We first study the 
finite range effects on the distorting potential potential 
in the incident channel following the method by John- 
son and Tandy [32] . In addition, we consider finite range 
effects in the evaluation of the transfer matrix element 
and the accuracy of the LEA. We also explore the impli- 
cations of our study to reactions involving loosely bound 
nuclei. 

The paper is organized in the following way: reaction 
theory is revised in Sec. II, results are presented in Sec. 
Ill and further discussed in section IV. Finally, conclu- 
sions are drawn in Sec. V. 



II. THEORY 

Our starting point is a three-body model of the n+p+A 
system, see 20], Ch.7 . In this model the scattering wave- 
function ^'•^•'(r, R) corresponding to a deuteron incident 
on a nucleus A is the solution of the inhomogeneous dif- 
ferential equation 

[E + ie-T,-Tn-UnA - UpA-Vnp]¥+\f,R) 

= ie(j)d{r)exp{iKd-R)- (1) 

Here T^ and Tr are the total kinetic energy operators 
associated with the n-p relative motion and the motion 
of the n + p centre of mass relative to the target, where 
r = Tp — r „ and R = (r„ + ^p)/2. We take r„ and r^ 
to be the neutron and proton coordinates relative to the 
center of mass of the target A. 

In Eq.([T|) the interactions UnA{rn), UpAirp) and Vnp{r) 
are the neutron-target, proton-target and neutron-proton 
interactions, respectively. The term proportional to le on 
the r.h.s. of Eq.([T|) ensures that there is an incoming wave 
only in the deuteron channel. 

One way of calculating the exact (d,p) transition am- 
plitude, T, from ^'^+^(r, ff^is to use the formulation of 
Timofeyuk and Johnson [2ll|. In the limit e — )■ 0, the 
transition amplitude T is given by: 



T=< 



.Axl~J\Vnp\^^+^>, 



(2) 



where (j)nA is the final state of the neutron-target system, 
and X g has a plane wave proton and an incoming scat- 
tered wave distorted by U*^ [2l|, [2j] . Note that in this 
formulation there is no remnant term and the proton dis- 
torted wave is generated by U*^ not U*g . Eq. ([4]) neglects 
recoil effects of order 1/A which were evaluated in [2ll |. 
They are negligible in all the cases discussed here except 



possibly for ^^C. The connection between this formula- 
tion and standard 3-body methods based on the Faddeev 
equations [2^ is explained in[23l|. 

A more common expression for the transition ampli- 
tude for this process is: 



T==< 



lAXpB 



v„. 



t/.n4-C/.Bl*^ + ^ 



JpA 



pB\ 



> 



(3) 



where (f)nA is the final state of the neutron-target system, 
and XpB i^ ^ proton scattering wave distorted by U*g. 
In the many-body generalisation of Eq.(j31) the remnant 
term UpA — UpB is a function of the internal coordinates 
of A and B and makes the interpretation of the tran- 
sition amplitude in terms of nuclear structure (overlap 
functions) much more complicated than when Eq.(l2]) is 
used. This introduces into the formulation an additional 
optical potential UpB and thus larger uncertainties into 
the analysis since in most applications to exotic nuclei 
this potential is not well determined. Many of the recent 
applications have used Eq. ([3]) and neglected the remnant 
term UpA - UpB^■ 



T=<0„AxyiKp|*(+)> 



(4) 



Neglecting the remnant term is a very good approxima- 
tion for all cases discussed here with the exception of ^^C. 
Eq.Q is the starting point for the present study. 

The important realization in 28, 32,] is that with Eq.® 
or Eq.Q, the full three-body wavefunction 'J'*^+^(r, ^) is 
only required within the range of the Vnp interaction. A 
major simplification is achieved in the limit of the zero- 
range approximation as then only ^'^+^(0,-/?) is needed. 



where the effective potential Ujs does not describe 
deuteron elastic scattering, but rather incorporates 
deuteron breakup effects within the range of Vnp'. 

Ujs{R) = UnA{R) + UpA{R). (8) 

Within this model, the transfer amplitude reduces to 

Tjs = Do J dR4>:,A{R)x;B{R)xf{R) (9) 

where Dq is the zero range constant of the deuteron. 



B. The Johnson and Tandy generalisation of the 
Johnson Soper method 

The Johnson and Tandy [S^] approach again builds on 
the fact that the three-body wavefunction is only needed 
within the range of Vnp- While the continuum discretized 
coupled channel method (CDCC) JS^ uses a basis of 
eigenstates of the n — p Hamiltonian which is complete 
for all values of the n — p separation f, as in Eq.([6]), here 
the Weinberg basis is introduced: 



{Tr + aiVnp)(j)i{r) = ~Si4>i{r) 



(10) 



with 



— 1,2,... and the orthonormality relation 
{4'i\Vnp\(j)j) ~ ^^ij- The Weinberg states (or Sturmi- 
ans) form a complete basis within the range of the Vnp 
interaction and thus they are particularly suited to de- 
scribing the problem when using the transfer amplitude 
written as in Eq. (j4]) . A clear advantage of this basis as 
compared to Eq.® is that it is square integrable. 
The three-body wavefunction is then expanded as: 



A. Johnson and Soper method 

The method developed by Johnson and Soper [28,] is 
based on an expansion of the three-body wavefunction 
^^"•"^ (r, R) in the complete set of eigenstates of the n—p 
Hamiltonian: 



{Tr + Vnp)4id{'r) = -ed(t>d{'r) , 
{Tr + Vnp)(t)k{r) = +ekMr)- 



(5) 



Here 4>d{r) is the deuteron wave function while scattering 
states are represented by 4>k{r). The three-body wave- 
function is then expanded as: 

¥+\r,R) = Mr)Xd{R) + f dk4^\r)xk{R) ■ (6) 



When this expansion is introduced in Eq.(IT]), and as- 
suming that the excitation energies of the deuteron are 
small compared to the deuteron initial kinetic energy 
Sk + Sd « E, the three-body equation for r = re- 
duces to an optical model equation: 



¥+\r,R)=J2M^MR)- 



fll) 



j=i 



{E + ed-TR- Ujs{R))xriR) = , 



(7) 



When this expansion is introduced into the three-body 
equation Eq.(IT|), one obtains a set of coupled channel 
equations: 

[E + ie - Kr - UuiRMiR) = 

ie6,iNd e^iKd-R) + Y. U^J{R)xAR)■ (12) 

The coupling potentials are defined by Uij{R) — Uij + 

Pijiaj - 1) and Uij{R) = -{(t)i\VnpiUnA + UpA\(t)j) where 
fHij = {4'i\V^p\(j)j) and aj are the eigenvalues of the Wien- 
berg equation Eq. (jlOp . The normalization coefiicient ap- 
pearing on the r.h.s of Ea. (IT2)) is Nd = — {(j)i\Vnp\(j)d} ■ 

These coupled channel equations can be solved exactly 
as done in 33] but reduce to a much simpler form if only 
the first term of the expansion Eg. ljlip is necessary. In 
that case, ai = 1 and we can arrive at the following 
optical model type equation: 

{E + Ed + ie-TB- Un{Rj)xi^{R) = *e^rf exp{iKd.R), 

(13) 



where now the potential is still related to the sum of the 
proton and neutron potentials as in Johnson and Soper, 
but involves a more complex folding procedure: 

UniR) = ~{Mr)\Vnp{UnA + UpA)\M^) ■ (14) 

Apart from the normalization, 0i is the ground state 
wavefunction of the deuteron cjjii. Then the potential 
Uii{R) can also be written in terms of (j)d'. 



{M^\Vnp(UnA + UpA)\M^) 

{Mr)\Vnp\Mr1) 



Uii{R) = Ujt{R) 

WdVllVnvl'PdVI'l/ 

(15) 

In this case the transfer amplitude is defined through 
the 6-dimensional integral: 



T = <<l)nAX^fJ\Vnp\(blx(^{R)> 

= < ^nAxi^^\Vnp\MxnR)/Nd) > 



(16) 



where we have used \ (j)d >= Nd\ (pi >■ 

We see from Eq.([T31) that {xi^{R)/Nd) is a distorted 
wave generated by the potential Un and normalized to 
an incident wave of unit amplitude. 

In the following section we will compare the cross sec- 
tions obtained with Eq.® and Ea.([TB)). We also dis- 
entangle the separate effects of finite range in the poten- 
tial f/ii, looking specifically at the potentials from Eq.® 
and ([T?l) , and that of finite range in the evaluation of the 
transfer amplitude. 



III. RESULTS 

We perform a systematic study of finite range effects 
on 26 (d,p) reactions. In all our calculations we take the 
Reid interaction for the deuteron ^] and the Chapel Hill 
glob al parameterization for the nucleon optical potentials 
|35| . In calculating the potential of EqUJlwe neglect the 
d-wave part of the (jyd- For all cases here presented, the 
final bound single particle state is obtained using a poten- 
tial with standard radius and diffuseness r = 1.2 fm and 
a = 0.65 fm and adjusting the depth to the known neu- 
tron separation energy of the corresponding final state. 
In Subsection(IIIA) we show our results for the potential 
and then present the results for the (d,p) cross sections 
in Subsection(IIIB). 



that the most important difference between the interac- 
tions UjT and Ujs is a constant increase in diffuseness. 
There is also a slight systematic decrease in radius. Dif- 
ferences in the depths of the real and imaginary parts are 



more subtle and vary case to case. In 37[, an approxi- 
mate method to estimate finite range corrections to the 
adiabatic potential was developed. In that method, the 
radius is fixed but an increase in diffuseness is predicted 
with a decrease in the depth of the potential. In table 
m we show the percentage difference of our numerically 
calculated Ujt and the Wales and Johnson approximate 
prescription Uwj [37|, relative to the zero-range Johnson 
and Soper potential Ujs^ for three reference cases. Per- 
centage differences are calculated at the radius for which 
the potential is maximum. 

The main feature of Ujt compared to Ujs is captured 
by the Wales and Johnson prescription, namely the in- 
crease in the diffuseness in both the real and imaginary 
part of the interaction. However the Wales and Johnson 
results differ quantitatively from ours. 



TABLE I: Finite range effects on the Johnson- Tandy distort- 
ing potential, Ea. (ll5|l : in the 3rd column the Wales and John- 
son potential is compared with the zero range potential Ujs 
and in the 4th column the Johnson and Tandy Ujt is com- 
pared with Ujs- We compare the diffuseness of the real and 
imaginary parts, an and a/, as well as the depths of the real 
and imaginary parts V and Ws, and the corresponding radii 
tr and r/. 



target 


parameter 


Uwj 


Utj 


all 


aR 


+4% 


+7% 




ai 


+3% 


+8-9% 


12c 


V 


-5.6% 


-1.98% 




TR 


0% 


-1.25% 




Ws 


-4.6% 


-4.52% 




ri 


0% 


+0.72% 


«Ca 


V 


-2.1% 


-0.04% 




TR 


0% 


-0.93% 




Ws 


-3.7% 


+1.6% 




ri 


0% 


-0.97% 


208pb 


V 


-0.7% 


+0.06% 




TR 


0% 


-0.35% 




Ws 


-3.3% 


+ 1.2% 




ri 


0% 


-0.35% 



A. Finite range effects in the potentials 



Johnson and Tandy potentials Eq. p5|) are computed, 
using a subroutine contained in the code TWOFNR [3a| 
for performing the r integrations needed, and compared 
with the Johnson and Soper adiabatic potentials Eq.® 
for 26 cases. In order to simplify the comparison we fit 
the real part of the resulting Ujt to a volume Woods 
Saxon shape and the corresponding imaginary to a sur- 
face Woods Saxon form. For all cases studied we find 



B. Finite range effects in the transfer cross sections 

Once the adiabatic potentials are defined, cross sec- 
tions can be obtained. The matrix element in Eq.® was 
evaluated using the adiabatic wavefunction distorted by 
Ujs and the zero-range constant Dq was obtained from 
the Reid n-p interaction [3^] for consistency. The fi- 



TABLE II: Percentage differences of finite range effects in 
(d,p) reactions relative to the zero-range Johnson and Soper 
calculation (* denotes cases for which no data is available). 
The target nucleus, the deuteron incident energy in the lab- 
oratory (in MeV) and the angle (in degrees) at which the 
percentage differences are evaluated are given in the 1st, 2nd 
and 3rd column respectively. 



Target 


Sd 


e 


A(LEA) 


A(FR-JS) 


A(FR-JT) 


A(JT-JS) 


12^ 


4 


25 


+5.6% 


+5.5% 


+4.5% 


-1.0% 


12^ 


12 


13 


+2.6% 


+2.9% 


-1.5% 


-4.3% 


12^ 


19.6 


10 


+11% 


+13% 


+7.7% 


-4.2% 


12^ 


56 


6 


-37% 


-27% 


-36% 


-12% 


«Ca 


2 


180 


+6.5% 


+6.3% 


+2.6% 


-3.5% 


«Ca 


13 


12 


+4.9% 


+3.8% 


-2.8% 


-6.2% 


«Ca 


19 


8 


+5.0% 


+4.0% 


-0.30% 


-4.1% 


«Ca* 


30 


4 




+7.3% 


+4.8% 


-2.3% 


«Ca* 


40 







-5.4% 


-5.9% 


-10% 


«Ca* 


50 







-1.9% 


-19% 


-18% 


«Ca 


56 





-5.2% 


-6.5% 


-24% 


-18.6% 


«9Ga 


12 


14 


+4.3% 


+4.7% 


-1.1% 


-5.49% 


«6Kr 


11 


25 


+4.8% 


+5.5% 


-0.40% 


-5.63% 


9°Zr 


2.7 


138 


+6.2% 


+7.3% 


+5.5% 


-1.7% 


3«Zr 


11 


26 


+5.4% 


+5.0% 


-0.90% 


-5.6% 


^^*Sn 


5.6 


175 


+6.1% 


+11% 


+7.5% 


-2.8% 


^^*Sn 


33.3 





+2.9% 


+4.6% 


0% 


-4.4% 


124gj^* 


40 


12 




-1.1% 


-2.4% 


-1.4% 


124s^* 


50 


11 




-3.9% 


-4.3% 


-0.44% 


124SJJ* 


60 


9 




-11% 


-30% 


-21% 


124gj^* 


70 





+5.1% 


-29% 


-44% 


-21% 


208p^^ 


8 


180 


+6.1% 


+7.2% 


+6.1% 


-0.96% 


208p^j 


12 


98 


+5.7% 


+8.8% 


+2.2% 


-6.1% 


208p^!|! 


20 


30 




+4.5% 


-2.3% 


-6.6% 


208pj^!|! 


40 


9 




+1.4% 


-6.9% 


-8.1% 


208p^=|= 


60 







+0.14% 


-8.8% 


-9.0% 


208p^!|! 


80 







-62% 


-86% 


-63% 



nitc range calculation Eq. ([161) ^^es Ujt for the adiabatic 
wavefunction. 

To assess the relative importance of the finite range 
effect in the adiabatic potential and that on the eval- 
uation of the transfer amplitude, we also perform a fi- 
nite range calculation using Ujs- Finally, given that the 
LEA method [29] has been widely used in the past, we 
also perform a calculation where Eq.([S]) is evaluated us- 
ing the Johnson and Soper adiabatic potential Ujs but 
making the local energy correction |29|. A ll calculations 
are performed using the code FRESCO [39j |. 

In tableiniwe quantify these effects for all the reactions 
studied. All percentage differences are calculated at the 
first peak of the angular distribution (with the exception 
of the sub-Coulomb examples for which percentage differ- 
ences are calculated at backward angles) and are relative 
to the Johnson and Soper approach Eq.®: 



• A (LEA) shows the effect of finite range in the eval- 
uation of the T-matrix when the LEA is used in 
conjunction with the Johnson and Soper model; 

• A(FR-JS) shows the effect of finite range in the 
evaluation of the T-matrix when fixing the adia- 
batic potential to Ujs', 

• A(FR-JT) shows the full finite range effects when 
finite range is taken into account properly both in 
the evaluation of the T-matrix and the adiabatic 
potential (Ujt)', 

• A(JS-JT) is the percentage difference between fi- 
nite range calculations using Ujt and Ujs and 
therefore shows the effect of including finite range 
in the adiabatic potential. 

In the table |ll] we indicate the angle at which the per- 
centage differences of cross sections were calculated. 

In addition to the full table, we also select some an- 
gular distributions that illustrate the various trends we 
observed, presented in FigsHJSl Each plot contains four 
lines: a dotted line corresponding to the zero-range John- 
son and Soper calculation, a dashed line corresponding 
to the local energy correction to the T-matrix calculation 
with the Johnson and Soper potential, the long-dashed 
line corresponding to a full finite range T-matrix calcu- 
lation where Ujs is used, and the full line is a finite 
range T-matrix calculation with the finite range adia- 
batic potential Ujt- Data is also presented whenever 
available, but only as an indication that the ingredients 
of our model are realistic and therefore the magnitude of 
the finite range effects reliable. It is not the purpose of 
this work to extract spectroscopic information for these 
systems. 

We first look at sub-Coulomb transfer reactions, which 
are usually rather insensitive to the nuclear optical po- 
tential. In FigUfa) we show 48Ca(d,p)49Ca at Eiab = 2 
MeV and in FigH^b) 208pb(d,p)209pb at Eiab = 8 MeV. 
In both cases, the effects of finite range are only a few 
percent (3% in ^^Ca and 6% in ^'^^Pb), and these are 
mostly due to the approximation in the evaluation of the 
T-matrix and not from the adiabatic potential. In this 
case, LEA is able to capture most of the finite range ef- 
fects. In the sub-Coulomb energy regime, for all cases 
studied, errors in using the Johnson and Soper potential 
with the local energy approximation are below 5%. 

Most of the available (d,p) data on stable systems 
was taken at energies above the Coulomb barrier for 
10-20 MeV deuterons. In Fig. [2Ja) we show results 
for 69Ga(d,p)^"Ga at Eiab = 12 MeV, in Fig. ^h) 
S6Kr(d,p)^'^Kr at Eiab = H MeV and in Fig. [^c) 
2"8pb(d,p)209pb at Eiab = 20 MeV. The overall effect 
of finite range in all three cases is very small (—1% for 
the 69Ga, 0.4% for the ^eRr and 6% in ^ospb), although 
it results from the cancellation of the two separate ef- 
fects, the finite range in the deuteron potential, which 
reduces the cross section and the finite range effect in 
the evaluation of the T-matrix which increases the cross 
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FIG. 1: Angular distributions for (d,p) at sub-Coulomb en- 
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and (b) ^°*Pb(d,p)^°^Pb(g.s.) £d = 8 MeV (data from ^). 
Comparison of zero-range Johnson and Soper model (dotted), 
the LEA Johnson and Soper model (dashed), a finite range 
calculation of the transfer amplitude using the Johnson and 
Soper adiabatic wave (long-dashed) and the full finite range 
results (solid line). 



section. No simple addition rule for these two effects was 
found. Here, the local energy approximation begins to 
show larger deviations from the full finite range calcula- 
tion. 

Finally, we also consider reactions at higher energies 
(50-80 MeV deuteron energy). Only two data sets are 
available, namely for ^^C and ^*Ca but we also include 
a study for ^^^Sn and ^°*Pb to ensure that our results 
are not biased by lower mass systems. All cases studied 
at these energies reveal that finite range effects are large 
and reduce the cross section. In Fig. [3lja) we show the 
angular distributions for ^^C(d,p)^'^C at Eiab = 56 MeV, 
in Fig. W!o) those for '*8Ca(d,p)49Ca at Eiab = 56 MeV 
and in Fig. [S^c) those for i24Sn(d,p)i24Sn at Eiab = 70 
MeV. The overall effect of finite range at the peak of the 
distribution is -36% for ^^C, -24% for "^Ca and -43.5% 
for ^■^^Sn. It is the finite range in the adiabatic poten- 
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tial that is the dominant cause for these large changes 
although the finite-range effect in the evaluation of the 
T-matrix is still important and should not be neglected. 
In addition, for these higher energies, we find that the 
LEA method breaks down and for the heavier systems 
this approximation can in fact provide a correction in 
the opposite direction to the full finite range calculation. 

Because it is the adiabatic scattering wavefunction that 
is mostly responsible for the large differences, we investi- 
gated the radial behavior of the scattering wavefunctions 
using either Ujt or Ujs for the partial waves which con- 
tribute the most to the transfer cross section. We specifi- 
cally looked at the properties of the integrand of Eq. p^ 
in the zero-range approximation where it has a simpler 
form. We found that the percentage difference comes 
from subtle cancellations and cannot be well illustrated in 
the partial wave expansion. Intuitively one might argue 
that since the energies are larger, the dominant contri- 
bution to the transfer cross section comes from smaller 
impact parameters and thus sensitivity to the range of 
Vnp should be larger. 

To ensure that our results are general, in particu- 
lar that they will still be applicable to reactions in 
which the final bound state has a large spatial exten- 
sion, we performed additional calculations for a fictitious 
^^Ca(d,p)''^Ca setting the valence neutron angular mo- 
mentum in the final bound state to ^ = state, and vary- 
ing the binding energy Sn — 0.1 — 6 MeV. The overall 
findings did not change: regardless of the loosely bound 
nature of the final nucleus, or the angular momentum in 
the final bound state, the effects of finite range in the 
transfer cross section are modest for low energies and 
become very important for the higher energies. 



IV. DISCUSSION 

The overall features obtained in this work can be best 
summarized in FigU] where the two separate effects of fi- 
nite range are plotted as a function of beam energy for 
(d,p) reaction on four different targets: solid symbols pro- 
vide the percentage effect of including finite range in the 
evaluation of the matrix element relative to a zero-range 
T-matrix calculation with a Johnson and Soper potential 
in both cases, and the open symbols correspond to the 
effect of including finite range effects in the distorting 
potential in the incident channel wave. The figure sum- 
marizes the results already given in Table |TT1 From this 
figure we can see that at low energy, finite range results 
differ by less than 10% from zero range matrix element 
with a Johnson and Soper adiabatic potential. However, 
as the incoming deutcron energy increases, finite range ef- 
fects become very important and can dominate the result. 
The energy at which the transition occurs depends non 
linearly on the charge and the mass of the system. For 
practical purposes we find the transition to be around 20 
MeV/u for lighter systems and 30 MeV/u for the heavy 
systems. 
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FIG. 3: Angular distributions for (d,p) at high energies: 
(a) i2C(d,p)"C(g.s.) Ed = 56 MeV (data from [H), (b) 
''*Ca(d,p)*^Ca(g.s.) Ea = 56 MeV (data from ^) and (c) 
^^''Sn(d,p)^25gjj(gg -) ^^ ^ 7Q ]y[gy Comparison of zero- 
range plus Johnson and Soper method (dotted), the LEA plus 
Johnson and Soper method (dashed), a finite range calcula- 
tion of the transfer amplitude using the Johnson and Soper 
adiabatic potential (long-dashed) with the full finite range 
results (solid line). 
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FIG. 4: Systematic finite range effect as a function of beam 
energy: open symbols give tlie effect in tlie incident cliannel 
distorted wave, and the filled symbols are the effect in the 
evaluation of the matrix elements. 



To facilitate the practical analysis of experiments, we 
searched for a global correction factor, a factor that 
would estimate the finite range effect as a function of tar- 
get charge, mass and beam energy. Although for a given 
target one could always find a function F{Ed) represent- 
ing the finite-range correction, no consistent dependence 
in mass and charge was found for the parameters of the 
various fits. 

It is important to remember that the exact inclusion 
of deuteron finite range effects requires the solution of a 
couple channel equations Eq. p^ and here we truncated 
the Weinberg expansion to the first term to simplify the 
problem. The full equations were solved in |33|] for ^^Zn. 
It would be interesting to extend this study to better 
determine the range of validity of the truncation here 
used. 



V. CONCLUSION 

We perform a systematic study of the effects of 
deuteron finite range in (d,p) reactions, within a formal- 



ism that includes the coherent effects of deuteron breakup 
through an adiabatic potential in the incident channel. 
We use the adiabatic formalism developed by Johnson 
and Soper [28] in zero range and compare with the finite 
range generalization of Tandy and Johnson |32] . We ana- 
lyze separately the effects of finite-range in the adiabatic 
distorting potential and finite range in the evaluation of 
the transfer matrix element. We also test the local en- 
ergy approximation which is widely used as an estimate 
of finite range corrections to zero range transfer cross 
sections. We performed (d,p) calculations to determine 
angular distributions for a wide range of beam energies 
as well as a variety of targets, from A=12 to A=208. 



For sub-Coulomb reactions, the percentage difference 
between the finite range and the zero range cross sections 
at the peak of the angular distribution relative to the 
zero range Johnson and Soper prediction is within 10% 
for all cases studied, and the local energy approxima- 
tion provides an estimate within a few percent of the full 
finite range calculation. However, as the beam energy in- 
creases, finite range effects become more important. For 
intermediate energies {E < 20 MeV/u for A < 50 and 
E < 30 MeV/u for heavier nuclei), including the finite 
range of the n — p interaction in the adiabatic scattered 
wavefunction reduces the cross section while including 
finite range in the evaluation of the transfer amplitude 
increases the cross section. Both effects are significant, 
although strong cancellations may occur. At higher en- 
ergies, both finite range effects have the same sign, re- 
ducing the transfer cross section. In this case we find the 
total effect of finite range to be very important. Our re- 
sults also suggest that at these higher energies, the local 
energy approximation is no longer adequate. 
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